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ABSTRACT 

The distance to the Large Magellanic Cloud (LMC) is crucial for the calibration of the Cosmic Distance 
Scale. We derive a distance to the LMC based on an analysis of ground-based photometry and HST-based spec- 
troscopy and spectrophotometry of the LMC eclipsing binary system HV2274. Analysis of the optical light 
curve and 7/5T/GHRS radial velocity curve provides the masses and radii of the binary components. Analy- 
sis of the HST/FOS UV/optical spectrophotometry provides the temperatures of the component stars and the 
interstellar extinction of the system. When combined, these data yield a distance to the binary system. After 
correcting for the location of HV2274 with respect to the center of the LMC, we find c/lmc = 45.7 ± 1.6 kpc or 
(Vq -M,,)lmc = 18.30 ± 0.07 mag. This result, which is immune to the metallicity-induced zero point uncertainties 
that have plagued other techniques, lends strong support to the "short" LMC distance scale as derived from a 
number of independent methods. 

Subject headings: Binaries: Eclipsing - Stars: Distances - Stars: Fundamental Parameters - Stars: Individual 
(HV2274) - Galaxies: Magellanic Clouds - Cosmology: Distance Scale 



1. INTRODUCTION 

We present the first accurate distance determination to the 
Large Magellanic Cloud (LMC) using an eclipsing binary sys- 
tem, resulting from an international collaboration studying the 
physical properties of these important objects. This program 
includes both ground-based and HST observations and is aimed 
primarily at determining the stellar properties (masses, radii, 
luminosities, and ages) and the distances of about a dozen se- 
lected +14 < V < +16 mag eclipsing binaries in the Clouds (see 
Guinan et al. 1996). It is well known that studies of eclipsing 
binaries yield the most fundamental determinations of the ba- 
sic stellar properties. We demonstrate here their great value as 
"standard candles" for distance determination. 

The distance to the LMC is a critical rung on the cosmic 
distance ladder and numerous independent methods (utilizing, 
for example, RR Lyrae stars, Cepheids, SN 1987a, and "red 
clump" stars) have been employed to determine it. Unfortu- 
nately, and despite the use of the eagerly-awaited Hipparcos 
parallaxes, there remains considerable disagreement about this 
important measurement and, as summarized by Cole (1998; see 
also Westerlund 1997), the current uncertainty is about 10-15%. 
In this paper, we show that the analysis of several eclipsing 
LMC binaries could reduce the uncertainty in the LMC distance 
to level of several percent. 

These first results deal with the LMC eclipsing binary Har- 
vard Variable 2274 (hereafter HV2274; < V > max ~ 14.2; Bl-2 
IV-III + Bl-2 rV-m; P = 5.73 days). This system is relatively 
bright, lightly reddened, and has complete CCD light curves 



showing it to have deep eclipses and uncomplicated out-of- 
eclipse light variations (Watson et al. 1992). Moreover, it is a 
detached system with nearly identical bright stars located in an 
uncrowded field. To provide an accurate distance to the system 
and, hence, the LMC, absolute radii, temperatures, and redden- 
ing corrections are needed (Guinan 1993; Gimenez et al. 1994; 
Guinan et al. 1996). The radii are obtained from the light and 
radial velocity curves and the temperatures and reddening from 
UV/optical spectrophotometry. The following sections present 
the HST and ground-based observations used in this study (§2), 
describe our analysis techniques (§3), and present our distance 
determination for the LMC along with a discussion of the un- 
certainties in this result (§4). 

2. OBSERVATIONS 

The ground-based V band light curve data for HV2274 are 
taken from Watson et al. (1992) and are shown in the lower 
panel of Figure 1 ("plus signs"). The light curve is well-defined 
and shows two minima of about equal depth (^0.7 mag) and 
nearly constant light outside of the eclipses. The more shallow 
minimum (secondary eclipse) occurs at orbital phase 0.53, in- 
dicating an eccentric orbit (<? = 0.136). We also adopt the very 
recently available V and B - V measurements from Udalski et 
al. (1998), i.e., V = 14.16 and B-V = -0.13. 

The spectrophotometry and radial velocity curves needed to 
complete the analysis of HV2274 were obtained by HST as part 
of our larger project on eclipsing binaries in the Clouds. Four 
FOS spectra, covering four different wavelength regions, were 
acquired, calibrated, and merged to produce a single spectrum 
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spanning 1150 A to 4820 A. These observations were made 
at orbital phases outside of the eclipses, when both stars were 
completely unobscured. 

We also obtained 16 HST/GERS medium resolution (R = 
23000) spectra at a number of selected orbital phases to 
construct the radial velocity curve. The spectra covered 
34 A and were centered at two different wavelengths: one near 
A1305 A and the other near A1335 A. These regions contain 
strong UV photospheric lines in early B stars (see Massa 1989). 
Stellar Fe III A1292.2, Si III AA1300 triplets (1294.5, 1296.7 and 
1298.9 A), and C II A1323.9 lines were fit with double gaus- 
sians, which gave the radial velocities for each component. The 
individual line measurements for each component were aver- 
aged to obtain its radial velocity at each phase. In the upper 
panel of Figure 1, these radial velocities are plotted against or- 
bital phase as determined by Watson et al. (1992; filled sym- 
bols). The mean uncertainties in the radial velocity measures 
are approximately ±15 km s" 1 . 

3. ANALYSIS 

3.1. Modeling the light and radial velocity curves 

The radial velocity and light curve data were analyzed using 
an improved Wilson-Devinney light curve analysis code (Wil- 
son & Devinney 1971; Wilson 1990) that includes Kurucz at- 
mosphere models for the computation of the stellar radiative 
parameters (Milone et al. 1994). The Wilson-Devinney code is 
a standard tool in the analysis of eclipsing binaries. 

The radial velocity and the light curve solutions for HV2274 
are shown in the top and bottom panels of Figure 1 , respectively 
(solid curves). For the light curve, the differences between ob- 
served and computed values at V ("(O-C)y") are also shown. 
Tables 1 lists the values determined for the orbital period (P), 
eccentricity (e), orbital inclination (;), longitude of periastron 
(uj), velocity semi-amplitudes (K), systemic velocity (7), tem- 
perature ratio (T B /T A ), luminosity ratio (L B /L A ), absolute stellar 
radii (r), masses (M), and surface gravities (g). 

The accurate measurement of the radii (r) of the components 
is critically important for deriving the distance to HV2274. The 
radii are determined by combining the fractional radii (rf = r/a) 
obtained from the light curve analysis with the orbital semi- 
major axis (a) derived chiefly from the spectroscopic orbit. 
Thus, the absolute radius of each star is r = rj x a and the cal- 
culated uncertainties in the radii are about 2.3%. The precision 
achievable in the measurement of fundamental properties, such 
as radius, is one of the keys to the usefulness of eclipsing bina- 
ries as standard candles. 

3.2. Modeling the spectrophotometry 

The observed energy distribution of HV2274, as ob- 
tained by the FOS between 1150 A and 4800 A is shown in 
Figure 2 (filled circles; the lower of the two full spectra). The 
observed fluxes depend both on the surface fluxes of the bi- 
nary's components and on the attenuating effects of distance 
and interstellar extinction. This relationship can be expressed 
as: 

/ Affi = l[r^ + 4^]xl0-°- 4 ^ (1) 
or, substituting for the total extinction, A(A), 

fx®={^)\F£ + (r B /r A ) 2 F*} (2) 
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where F' x {i = A,B} are the emergent fluxes at the surfaces of the 
two components of the binary, the r, are the radii of the com- 
ponents, and d is the distance to the binary. A(A) is the total 
extinction along the line of sight to the system at each wave- 
length A and is expressed in Eq. 2 as a function of E(B - V), 
the normalized extinction curve k(X-V) = E(X-V)/E(B-V), 
and the ratio of selective to total extinction in the V band 
R = A(V)/E(B — V). Our analysis consists of modeling the 
observed energy distribution of HV2274 via a non-linear least 
squares procedure to derive values of the parameters (r A /d) 2 , 
F{, E(B-V), and k(X-V) in Eq. 2. Ultimately, the distance to 
the binary will be determined from (r A /d) 2 . 

We represent the stellar surface fluxes F l x by ATLAS9 atmo- 
sphere models from R.L. Kurucz, which each depend on four 
parameters: effective temperature (T e ff), surface gravity (g), 
metallicity (normally z), and microturbulence velocity {p). The 
eight parameters needed to define the surface fluxes of the two 
binary components are constrained by the results of the light 
and radial velocity curve analyses and we adopt the effective 
temperature ratio T B /T A and \ogg values listed in Table 1 . In ad- 
dition, we assume that both components of the binary have the 
same metallicity and the same microturbulence velocity. The 
ratio of the stellar radii (r B /r A ) 2 in Eq. 2 is determined from the 
light and radial velocity curve solutions (see Table 1), and we 
adopt the standard Galactic value of R = 3.1 (see §4). 

The normalized extinction curve k(X - V) is modeled using 
the results of Fitzpatrick & Massa (1990), who showed that the 
functional form of the extinction wavelength dependence in the 
UV is tightly constrained. The parameters describing the shape 
of the curve are determined from the fitting procedure. The 
method of smoothly joining UV and optical extinction curves 
is discussed by Fitzpatrick (1999). 

This technique of solving simultaneously for both the stellar 
parameters of a reddened star and the shape of the UV/optical 
extinction curve is described in detail by Fitzpatrick & Massa 
(1999; FM99). They demonstrate that the new ATLAS9 models 
can reproduce the observed UV/optical continua of unreddened 
main sequence (and slightly evolved) B stars to a level consis- 
tent with the uncertainties in currently available spectrophoto- 
metric data. They further show that both the model atmosphere 
parameters and the wavelength dependence of interstellar ex- 
tinction can be extracted from analysis of UV/optical spectra of 
reddened stars. This is possible because the spectral "signature" 
of interstellar extinction is very different from the "signatures" 
of temperature, surface gravity, metallicity, or microturbulence. 

The dereddened energy distribution of HV2274 is shown as 
the upper spectrum in Figure 2 (small filled circles with obser- 
vational error bars) superimposed with the best-fitting model 
(solid "histogram" curve). The inset shows a blowup of the 
Balmer jump region. The dereddened V and B data from Udal- 
ski et al. (1998), converted to flux units, are shown by the 
large filled circles. The excellence of this fit (the reduced x 2 
is close to 1) is illustrative of the fits achieved by FM99 for B 
stars in general. The parameters of the fit are listed in Table 1. 
For B-type stars, the "metallicity" measured from modeling the 
UV/optical flux is most heavily influenced by absorption from 
iron-group elements in the UV region. Therefore in Table 1 we 
refer to the derived metallicity as [Fe /H], and the result is quite 
reasonable for a LMC star. The uncertainties indicated in Ta- 
ble 1 are \-a internal errors and incorporate the full co variance 
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(or interdependence) of all the parameters. I.e, if any of the 
parameters is changed by ± 1-cr, the total x 2 of the best fit — 
after all the other parameters are reoptimized - increases by + 1 . 
The small size of these uncertainties is testament to the lack of 
covariance among the parameters and the quality of the data. 

The best-fit model found in this paper differs significantly 
from that reported in an earlier version of this work (Guinan et 
al. 1998). This results entirely from the inclusion here of the V 
and B — V measurements from Udalski et al. (1998). The ear- 
lier work incorporated very uncertain optical photometry from 
Watson et al. (1992), but these data were weighted so low that 
results were based essentially entirely on the FOS data. This 
had an important effect because, unless the observations extend 
through the 4400 - 5500 A region, a possible degeneracy ex- 
ists between the best-fit values of E(B - V) and the normalized 
extinction curve k{ A- V). I.e. very similar relative extinction 
corrections, E(B-V) x k(X-V), can be obtained with a range 
of E(B-V) values. In Guinan et al. (1998) the degeneracy re- 
gion was found to be E(B-V) ~ 0.08-0.12, with the best fit 
given by E(B-V) ~ 0.08 — the solution reported in that pa- 
per. This result appeared reasonable since, combined with the 
Watson et al. photometry (B-V = -0.18), it yielded an intrin- 
sic color of (B-V)o ~ -0.26, compatible with the spectral type 
of HV2274. Using the more accurate (and thus more heavily 
weighted) Udalski et al. photometry, the reddening degeneracy 
disappears and a reasonable fit to the entire observed spectrum 
can only be achieved with E(B-V) = 0.12 — a value at the 
opposite end of the degeneracy "valley" from that reported in 
Guinan et al. (1998), as expected from Murphy's Law. 

4. THE DISTANCE TO THE LMC 

The distance to HV2274 is computed from the values of 
(r^/aT) 2 and r& derived above, and the uncertainty in this re- 
sult is found by propagating the uncertainties in each quantity. 
The total uncertainty in (r^/fif) 2 has internal and external com- 
ponents. The internal 1-cr fitting error is listed in Table 1. The 
external effects include uncertainty in the the extinction ratio 
R, which we take to be l-a(R) = ±0.3, and uncertainty in the 
zero point of the FOS spectrophotometry, which we assume to 
be l-a(FOS) = ±2-3% (Bohlin 1996; Bless & Percival 1996). 
Quadratically combining these independent influences yields a 
total uncertainty of a[(rA/d) 2 ] ~ 4.5%. The uncertainty in 
is taken from the Wilson-Devinney analysis and is listed in Ta- 
ble 1. From these results we derive o?hv2274 = 46.8 ±1.6 kpc, 
or, (Vo-M,,)hv2274 = 18.35 ± 0.07 mag. This result is 0.14 mag 
smaller than that reported by Guinan et al. (1998) due to the 
resolution of the aforementioned reddening degeneracy. 

To obtain the distance to the center of the LMC, the above re- 
sult must be corrected for the position of HV2274 in the LMC. 
The line-of-sight to the HV2274 binary system projects onto the 
end of the western region of the central bar of the LMC. The 



//5T/GHRS spectra show strong ISM lines from LMC gas in 
O I A 1302.2 A and Si II A 1304.4 A. These lines show only one 
component with a Heliocentric radial velocity of 285 km s" 1 , 
similar to the HV2274 systemic velocity of 312 km s . These 
data suggest that HV2274 lies in or close to the disk of the 
LMC. To transform the distance of HV2274 to the LMC opti- 
cal center we assumed the position angle of the line of nodes 
is 168° and that the disk inclination is 38° (Schmidt-Kaler 
& Gochermann 1992). Simple trigonometry then implies that 
HV2274 is located on the far-side of the LMC, about 1 100 pc 
behind the center. The distance of the center of the LMC is then 
d LMC = 45.7 ± 1.6 kpc or (V -M v )uac = 18.30 ±0.07 mag. 

This result lends strong support to the "short" LMC distance 
scale (see Cole 1998). Its importance, however, lies not only in 
the value of the distance itself, but in the precision that appears 
possible from using eclipsing binaries. This precision can be 
understood by noting two facts: 1) the values of \ogg for the 
two binary members are very well determined by the light and 
radial velocity curve analysis; and 2) the best fitting model at- 
mosphere reproduces the HV2274 energy distribution over the 
entire observed range (1200 - 5500 A) to within the small ob- 
servational uncertainties (see Fig. 2). These imply that the stel- 
lar temperature must — like the gravity — be well-determined, 
since features such as the Balmer jump are sensitive to both 
properties. Further, a correct temperature could not be found 
unless the distorting effects of interstellar extinction were prop- 
erly estimated and removed, since the possibility of a degener- 
acy in the reddening solution has been eliminated by the inclu- 
sion of V band data in the analysis (§3.2) . Finally, we note that 
this technique is immune to the metallicity-induced zero point 
uncertainties that have plagued many other LMC distance esti- 
mates because the determination of the stellar metallicity is an 
explicit part of the analysis and because the derived distance is 
actually extremely insensitive to the metallicity. 

The result from a single binary system — as precise as it ap- 
pears to be — does not by itself resolve the LMC distance issue. 
Unanticipated external systematic effects, such as a peculiar lo- 
cation of the star within the LMC, could compromise the de- 
rived distance. Such effects can only be addressed through the 
analysis of more systems, preferably with a variety of locations 
within the LMC and covering a range of stellar properties. The 
ultimate precision of the LMC distance will best be evaluated 
by the range of results derived from such analyses. The results 
presented here, however, demonstrate the power of eclipsing bi- 
naries to address fundamental astrophysical issues, such as the 
derivation of basic stellar properties, and cosmologically im- 
portant issues, such as the LMC distance. 

This work was supported by NASA grants NAG5-7113 and 
HST GO-06683, and NSF/RUI AST-93 15365. 
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Table 1 

Orbital and Stellar Parameters 



Orbital Solution and Stellar Properties' 1 

P= 5.726006(12) days e = 0.136(12) 
; = 89°.6(1.3) a = 38.58(93) R Q 

w(1990.6) = 73°.3(1.5) 7 = +312(4) km s" 1 
K A = 166.2(5.9) km s" 1 K B = 177.3(5.8) km s" 1 
T B /T A = 1 .005(5) L B /L A (V) = 0.844(5) 

r A = 9.84(24)fl Q r B = 9.03(20)fl Q 

M A = 12.1(4)M Q M B = 11.4(4)M 

logg A = 3.54(3) logg B = 3.58(3) 

Stellar Properties from UV/Optical Spectrophotometry 

T A = 23000(180)/: T B = 231 10(1 S0)K 

[Fe/H] AB = -0.45(6) ^ab = 1 -9(7) km s" 1 
E(B-V) = 0.120(9) mag (r A /d) 2 = 2.249(63) x 10" 23 
Distance to HV2274 

d = 46.8(1.6) kpc Vp-M v = 18.35(0.07)mag 

Distance to Center of LMC 

d = 45.7(1.6) kpc V -M v = 18.30(0.07) mag 



a The uncertainties in the parameters resulting from the 
light curve and radial velocity curve analyses were set at 
three times the r.m.s. scatter of several solutions obtained 
from different initial conditions. 

b The uncertainties in the parameters derived from anal- 
ysis of the UV/optical spectrophotometry are 1-er internal 
errors derived from a non-linear least squares fitting proce- 
dure (see §3.2). 
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FIG. 1. — The upper panel shows the radial velocity measurements of the two components of the HV2274 binary system as derived from medium resolution 
HST/GHRS spectra (filled symbols) and the radial velocity curve solution derived from the Wilson-Devinney analysis (solid curves). The discontinuous jumps seen 
in the model curves, known as the "Rossiter effect," occur when rotating stars are partially eclipsed. The lower panel shows the ground-based light curve data from 
Watson et al. 1992, based on differential photometry (Variable - Comparison) in the V band ("V-Cp"; plus signs). The fit to the light curve is shown (solid curve) 
as well as the Observed - Computed residuals ("O-C"). 
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FIG. 2. — The UV/optical energy distribution of HV2274, in units of ergs cm~ 2 s~' A -1 . The lower full spectrum shows the observed HST/FOS energy distribution; 
the upper spectrum shows the extinction-corrected energy distribution, superimposed with the best-fitting ATLAS9 atmosphere model (plotted in histogram style). 
Vertical lines through the data points indicate the l-cr observational errors. Crosses indicate data points excluded from the fit due to contamination by interstellar 
absorption lines. The large filled circles show dereddened B and V photometry from Udalski et al. 1998. An expanded view of the fit near the Balmer jump is shown 
within the inset. 



